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PROCEDE DE RECONSTRUCTION D'UNE IMAGE RADIOGRAPHIQUE 
PAR COMBINAISON DE VIGNETTES SE RECOUVRANT 

DESCRIPTION 

5 

L' invention concerne un proc^de de 
reconstruction d'une image radiographigue par 
combinaison d'une collection de vignettes se 
recouvrant . 

10 Les images de grande taille en radiographie 

sont frequemment obtenues par morceaux, au moyen de 
vignettes de projection qui sont prises success ivement 
en d6pla<?ant le rayonnement traversant l'objet & 
diffSrentes positions, ainsi que le reseau de' 

15 detecteurs prenant les mesures. L' assemblage des 
vignettes .donne ensuite 1' image recherch^e. 

Le procede se complique dans le cas usuel 
d f un rayonnement divergeant d'un foyer vers, le reseau 
de detecteurs , soit en cone, soit en un ensemble 

20 d'eventails plans et parall^les. La figure 1 montre la 
configuration normale des mesures : le rayonnement 
comprend une source 1 (ponctuelle ou lineaire) qu f on 
deplace a chaque mesure le long de l'objet 2 ainsi que 
le reseau de detecteurs 3 „ Les positions prises sont 

25 notees la, lb, lc, Id, et 3a, 3b, 3c et 3d. Pour que 
1 1 attenuation du rayonnement soit mesur^e a tout 
endroit de l'objet 2, le faisceau 4 du rayonnement doit 
comprendre des portions de recouvrement dans les 
positions 4a, 4b, 4c et 4d qu'on lui fait prendre, et 

3 0 les vignettes de projection du rayonnement, dont les 
positions .coincident avec celles 3a, 3b, 3c et 3d que 

B 14215. 3/ JCI 
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le rlseau de detecteurs 3 prend successiveraent , 
presentent de meme des portions de recouvrement . II est 
done impossible de juxtaposer simplement les vignettes 
pour obtenir 1 1 image globale de l'objet, mais on doit 
5 au contraire determiner les positions des portions de 
recouvrement sur les vignettes et faire une synthese du 
contenu de ces portions de recouvrement pour 
reconstruire 1 1 image . 

Un autre probleme qui apparait est celui du 
10 grossissement des details selon leur distance a la 
source 1. La largeur de projection de details 5 de 
l'objet 2 sur le reseau de detecteurs 3 sera 
proportionnellement plus large si les details 5 sont 
plus proches de la source 1. Un rayonnement divergent 
15 ne permet ainsi ni de juxtaposer facilement des 
vignettes/ ni de respecter l'echelle des details a 
1 ' interieur de chaque vignette . 

La figure 2 permet d'expliciter ces 
problemes. Deux details 5a et 5b espaces verticalement 
20 se trouvent dans la portion de recouvrement des 
vignettes prises par le reseau de detecteur 3 aux 
positions 3a et 3b. Les rayons passant par le detail 5a 
sont distants de 1 ' ecart 6 sur le plan du reseau de 
detecteurs 3, et ceux qui traversent le detail 5b sont 
25 distants de 1 ' ecart 7 sur le meme plan ; les ecarts 
sont differents entre eux, et differents du deplacement 
qu'il a fallu imposer au reseau de detecteurs 3 entre 
les positions 3a et 3b ou les vues ont ete prises. Une 
bonne reconstruction de 1' image aux portions de 
30 recouvrement impose de combiner entre elles les mesures 
associees a chacun des details pour des vignettes 
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differentes, ce qui est impossible a faire directement 
puisque leurs hauteurs sont inconnues en general . Si on 
choisit par exemple d'associer les rayons distants de 
1'ecart 6 pour reconstruire les portions de 
5 recouvrement, les details a hauteur de 5a seront rendus 
correctement, mais les details presents a d'autres 
hauteurs ne pourront pas l'etre. La combinaison des 
vignettes produira alors du flou et un gross issement 
inexact de ces autres details. 

10 Un proc^dS de ce genre a cependant dej£ ete 

propose dans l'art. L' image est reconstruite en 
choisissant une reconstruction exacte des portions de 
recouvrement a des hauteurs d^terminees, la ou les 
details importants, et notamment les os pour vine-. 

15 radiographic du corps, sont susceptibles de se trouver. 
Pour obtenir des bons r€sultats, il faut done connaitre 
d'abord la hauteur de ces details. Des coefficients de . 
ponderation peuvent privilegier les resultats de l'une 
ou 1' autre des vignettes selon la position consideree 

20 sur la portion de recouvrement. La restitution des 
autres details de 1' image est sacrifice. 

Un autre precede connu consiste a calculer 
des correlations entre les portions de recouvrement des 
differentes vignettes pour evaluer 1'ecart (6, 7 ou 

25 autre) des rayons a associer pour faire la synthese des 
portions de recouvrement. Les correlations dependent de 
details preponderants presents sur les deux portions de 
recouvrement et provenant d'un meme endroit de 1'objet 
2. La reconstruction de 1' image est accomplie a hauteur 

30 de ces details preponderants et elle est bonne, si du 
moins ces details existent ; mais comme dans le procSde 
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precedent, les details situes aux autres hauteurs 
seront mal rendus, 

II faut a j outer que des conflits de hauteur 
peuvent apparaltre si les portions de recouvrement sont 
5 nombreuses, et en particulier avec un rayonnetuent 
conique ou les portions de recouvrement c one e men t tout 
le perim^tre des vignettes. Deux portions de 
recouvrement sur deux cotes d'une vignette pourront 
§tre reconstruites independamment a des hauteurs 
10 diff§rentes, tout en ayant une intersection pour 
laquelle on sera embarrasse a choisir une hauteur de 
reconstruction . 

Un procede de reconstruction plus correct e 
d'une image radiographique est propose avec 
15 1' invention. II est fonde- sur tine discretisation 
generale de 1'objet en volumes (voxels) definissant des 
hauteurs de reconstruction, et des corabinaisons des 
valeurs d' attenuation estimees sur chacun des volumes 
aux diffSrentes hauteurs de reconstruction pour 
20 am^liorer 1' image globale, sans privilSgier 
necessairement une hauteur de reconstruction. 

De fa<?on plus detaillee, 1' invention 
conceme generalement un proc€d£ de reconstruction 
d'une image radiographique d'un objet traverse par un 
25 rayonnetuent divergent subissant une attenuation, le 
rayonnement occupant des positions successives ayant 
des portions de chevauchement et 1 1 attenuation etant 
mesuree par ion reseau de detecteurs, sur lesquels le 
rayonnement se projette et dormant des vignettes de 
30 I 1 image respectivement associees aux positions du 
rayonnement et comprenant aussi des portions de 
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chevaucheraent , le proctde comprenant une combinaison de 
vignettes pour reconstruire 1' image, ainsi que les 
etapes suivantes : 

discrStiser l'objet en volumes 
5 definissant des hauteurs de reconstruction, 

associer les volumes a au moins un 
detecteur respectif du reseau sur lequel le rayonnement 
se projette apres avoir traversg ledit volume, 

attribuer une valeur d' attenuation a 

10 chaque volume d" apres les valeurs mesurees par ledit 
detecteur associe, 

- et combiner les valeurs d' attenuation des 
volumes aux differentes hauteurs de reconstruction. 

Dans une des formes de 1' invention, la 

15 valeur d' attenuation attribute a chaque volume est 
egale a la somme des valeurs mesurees par.' * ledit' 
detecteur associe, divisee par le nombre des vignettes 
que contribuent a donner ledit detecteur associe, et. 
les valeurs d' attenuation des volumes sont combinees 

20 par une combinaison numtrique sur des groupes des 
volumes superposes aux differentes hauteurs de 
reconstruction. Et dans une autre de ses formes, la 
valeur d' attenuation attribute a chaque volume est 
obtenue par retroprojection iterative des valeurs 

25 d' attenuation mesurees par les detecteurs, des valeurs 
provisoires etant attributes aux volumes et corrigees 
apres avoir ete projetees sur les ' detecteurs, en 
calculant des differences entre des sommes et les 
valeurs provisoires sur des lignes de projection aux 

30 valeurs mesurees par les detecteurs sur lesdites lignes 
de projection, et en repartissant les differences sur 



B 14215. 3/ JCI 



6 



lesdites lignes de projection pour corriger les valeurs 
provisoires. 

L' invention sera maintenant decrite 
complStement en liaison aux figures, dont la figure 1 
schematise le procede de realisation des vignettes, la 
figure 2 illustre le probleme de reconstruction a une 
hauteur arbitraire, la figure 3 illustre les elements 
participant a 1 ' exploitation de 1' invention, et les 
figures 4 et 5 sont des organigrammes de deux modes du 
proc€d€ . 

On passe au coramentaire de la figure 3. 
L'objet 2 est discretise en volumes ou voxels 8 qui 
definissent des hauteurs de reconstruction 11. Le 
rayonnement passe a travers les volumes 8 par des 
rayons 9i et 9 j , qui sont plusieurs pour les volumes 8 
appartenant aux portions de recouvrement , et qui sont 
originaires de positions respectives li et lj de la 
source 1 et se projettent sur des detecteurs respect if s 
lOi et lOj qui leur sont associ£s pour les positions 
correspondantes 3i et 3j du reseau de detecteurs 3. Les 
detecteurs 10 mesurent des attenuations des rayons 9 a 
travers tout l'objet 2, et done a travers tous les 
volumes tels que 8 qu' ils traversent. Dans la pratique, 
les volumes 8 se projettent sur une surface qui peut 
englober plusieurs detecteurs 10 completement , et 
d'autres partiellement * Le systeme est calibre pour 
associer a chaque volume 8 sa surface de projection et 
calculer l f attenuation qu' il fait subir au rayonnement 
en proposition des detecteurs concernes. Nous 
n'evoquerons pas ici ces techniques de calibration, qui 
sont tout a fait usuelles dans l'art, et considSrerons 
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des volumes 8 se projetant completement sur un 
detecteur 10 unique le long d'un rayon de projection 
unique par souci de simplicity des explications. 

D'apres la figure 4, le procede commence 
5 done par une 6tape A de discretisation generale de 
1'objet 2 en volumes 8 dont les couches definissent les 
hauteurs de reconstruction 11 de 1' image. En pratique, 
les hauteurs de reconstruction 11 seront assez peu 
nombreuses et les volumes 8 plut6t parallel^pipediques, 

10 allonges en hauteur, que cub i que s . Les etapes suivantes 
B et C consistent & se placer a une hauteur de 
reconstruction 11 et un volume 8 determines. On 
recherche alors les rayons tels que 9i et 9j passant.; 
par le volume 8 consid^re, et les detecteurs lOi et 10 j 

15 de projection desdits rayons, sur le reseau de 
detecteurs 3, a 1'etape D. L'etape suivante E consisted 
a lire la raesure d' attenuation des rayons 9i et 9j sur; 
les detecteurs lOi et 10 j . A l'6tape suivante F, une 
moyenne de ces attenuations est faite, du raoins pour-^ 

20 les volumes 8 appartenant aux portions de recouvrement 
et qui sont done traverses par au moins deux rayons 9. 
Les etapes C h E ou F sont ensuite repStees pour tous 
les volumes de la couche consideree ; apres quoi, a 
l'<§tape G, une image de 1'objet 2 est reconstruite . 

25 Cette image est une image de 1' ensemble de 

l'objet 2, et pas seulement une image de coupe a la 
hauteur consideree, puisque les attenuations mesurees 
par les detecteurs 10 le long des rayons 9 ont ete 
supposees concentrees aux volumes 8 de la couche & 

3 0 cette hauteur. 



B 14215. 3/JCI 



8 



Ensuite, on revient a 1'etape B pour 
reconstruire l'objet 2 a une autre hauteur, et le cycle 
des etapes C a G recommence avec les volumes 8 de la 
couche associee. 

Quand les images de 1 ' ob j et 2 ont 6t § 
reconstruites a toutes les hauteurs, elles sont 
combines & I'Stape H avec l'espoir d'obtenir une image 
plus exacte. Plusieurs procedes peuvent etre envisages* 
Le plus simple consiste peut-etre §l faire des moyennes 
des images sur des colonnes 12 {a la figure 3) de 
volumes 8 empilSs appartenant a des couches 
differentes, avec eventuellement une ponderation pour 
favoriser les couches les plus representatives, 
Eventuellement, on peut choisir une seule des images 
qu'on juge meilleure que les autres, ou un assemblage 
de plusieurs des images aux endroits qu' elles 
repr§sentent le mieux, Toutes ces methodes devraient 
donner de meilleurs rSsultats que celles de l'art 
anterieur qu'on a decrites auparavant . 

On ne fera que mentionner certains procedes 
de correction qui sont usuels dans l'art et qui ne sont 
pas affectes par 1' invention. 

Le rayonnement diffus6 peut tout d'abord 
etre soustrait des mesures avant d' exploiter celles-ci. 
Plusieurs procedes existent pour faire cette 
soustraction, dont le plus simple est peut -etre de 
proceder a une mesure supplementaire ou un Scran est 
intercalS entre l'objet 2 et le reseau de detecteurs 3 
en masquant certains des detecteurs 10- Les detecteurs 
10 masques ne sont pas touches par le rayonnement 
direct des rayons 9, mais seulement par le rayonnement 
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diffuse, qui est alors mesurg par ces detecteurs et qui 
peut §tre deduit par des interpolations pour les autres 
detecteurs . 

Les attenuations d'un rayonnetnent peuvent 
5 en general §tre exprimees par un coefficient 
multiplicateur du rayonnement initial Io infer ieur a 
1' unite et Sgal a e -|Lll, ou 1 est la longueur 
d' attenuation et \i le coefficient d' attenuation 
caracteristique du matSriau, et qui est en general la 

10 valeur qu'on recherche pour reconstruire 1' image, Les 
detecteurs 10 mesurent directement le rayonnetnent I qui 
n'a pas £te absorbe par l'objet 2 et qui est egal a 
Io e"^; on peut en deduire le produit p.1, puis la^ 
valeur de jx si on divise les valeurs de |ll par les 

15 longueurs de traversee de l'objet 2 par les rayons 9,' 
apr£s les avoir estim^es par une autre mesure ou les 
avoir evaluees. 

Un autre de mode de realisation de 
1' invention sera maintenant decrit au moyen de la' 

20 figure 5. Apr£s urie etape de discretisation J semblable 
a celle A de la realisation prgcedente, on effectue au 
mieux une division en blocs a 1' etape K. En effet, la 
resolution qui va etre entreprise peut devenir malaisee 
si le systeme considere est trop volumineu^ . En 

25 pratique chaque bloc pourra comprendre les volumes 8 
associes a une vignette. Qu'une division en bloc soit 
faite ou non, le probleme a resoudre peut s'exprimer 
par p = Mx ou x designe les inconnues, c'est-a-dire 
les attenuations aux volumes 8, p designe les 

30 projections de ces valeurs, c'est-a-dire les mesures 
par les detecteurs 10, et enfin M designe la matrice de 
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projection. Les coefficients mij de la mat rice M 
representent la contribution d'un volume 8 d' indice j a 
la projection suivant le rayon 9 d' indice i, et peuvent 
en general etre approches par la longueur traversEe par 
ce rayon dans ce volume . 

L'etape suivant e L est une evaluation de 
1' att§nuation aux volumes 8 du bloc considere. La 
premiere Evaluation peut dtre arbitraire, par exemple & 
des valeurs nulles * Pour chacun des volumes 8 on 
recherche, d'aprds l'Etape M, le detecteur 10 qui lui 
est associe par le rayon 9 qui le traverse, comtne & 
l'etape D de la realisation precedente . L'etape 
suivante N est une lecture des mesures des detecteurs 
10 semblable & l'etape E. La determination des rayons 9 
de projection permet de proceder k une evaluation des 
valeurs projetees d' attenuation a l'etape O, c'est-a- 
dire gu'on procede au calcul Mx pour ^valuer p. En 
soustrayant ces valeurs evaluees des projections aux 
valeurs reelles, raesurees, des memes projections, on 
determine l'erreur commise dans 1' Evaluation des 
valeurs projetEes a l'etape P. 

L'etape suivante Q est tone retroprojection 
de cette erreur dans les volumes 8 de l'objet 2 afin de 
corriger les valeurs evaluees de 1' attenuation. 
Concrltement , on procdde en executant la formule 

2 (q + D = 2 «D + X (q)|_!!W fe loc -M bloc ic ( *>) oil x<* +1) et x<* sont 

des evaluations successives de 1' attenuation aux 
volumes 8 du bloc ; X (q) est un coefficient de 
relaxation permettant de ne pas se diriger trop vite 
vers une solution qui ne correspond qu'aux premiers 
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blocs et qui est compris entre 0 et 2 ; ce coefficient 
n'est d'ailleurs pas uniforme dans les blocs mais peut 
avantageusement Stre plus el eve pour les rayons 
sensiblement verticaux, ou perpendiculaires aux 
5 detecteurs 10, afin de leur donner une importance plus 
grande ; ^bioc ^st la transposee de la matrice M pour le 
bloc considere ; le terme au denominateur est un terme 
de normalisation ; enfin, les termes entre parentheses 
representent l'erreur calculee a l'etape P . 

10 On procede de merae pour le bloc suivant, en 

recommengant le cycle de l'etape K & l'etape Q, puis on 
revient au premier bloc pour une nouvelle iteration, 
jusqu'a ce que les attenuations evaluees aient converge 
vers une solution, ce qu'on exprime par l'etape R. Les 

15 , volumes 8 inclus dans les portions de recouvrement du 
volume 2 ont ete traites de la mSrae fagon que , les 
autres, en subissant simplement des iterations plus 
nombreuses qu'eux s'ils sont associes S plusieurs blocs 
de projections, 

2 0 On dispose alors d'une image 

tridimensionnelle de l'objet 2 ; une image 
bidimensionnelle de bonne qualite peut §tre obtenue par 
une combinaison des valeurs obtenues, qui consiste a 
additionner les valeurs d' attenuation sur les colonnes 

25 12 de volumes 8 empiies. 

Le procede de 1' invention permet de 
concilier une bonne qualite de restitution des details 
importants de l'objet etudie avec une bonne qualite 
d' ensemble de 1' image, II est possible d'obtenir des 

30 images dont la resolution est analogue au pas des 
detecteurs 10 ♦ 
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On s'est place dans la situation usuelle ou 
le r^seau 3 de d€tecteurs accompagne le mouvement du 
rayonnement 4, mais le procede pourrait etre applique 
sans changement avec un reseau de detecteurs immobile 
sous 1'objet 2 et dont la surface s'etendrait a toutes 
les vignettes de projection. 
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REVINDICATIONS 

1) Precede de reconstruction d'une image 
radiographique d 1 un obj et traverse par un rayonnement 
5 divergent subissant une attenuation, le rayonnement 
occupant des positions successives (4) ayant des 
portions de chevauchement et 1 r attenuation etant 
mesur^e par un reseau (3) de detecteurs (10), sur 
lesquels le rayonnement se projette et dormant des 

10 vignettes de 1 ! image respectivement associ^es aux 
positions du rayonnement et comprenant aussi des 
portions de chevauchement, le procede comprenant une 
combinaison de vignettes pour reconstruire 1' image, 
ainsi que les Stapes suiyantes : 

15 - discretiser l'objet en volumes (8) 

def inissant des hauteurs de reconstruction (11) , 

associer les volumes a au moins tin 
detecteur respectif du reseau sur lequel le rayonnement 
se projette aprSs avoir traverse ledit volume, 

2 0 - attribuer une valeur d 1 attenuation a 

chaque volume d'apres les valeurs mesurees par ledit 
detecteur associe , 

- et combiner les valeurs d 1 attenuation des 
volumes aux diff ^rentes hauteurs de reconstruction, 

25 2) Proced6 de reconstruction d f une image 

radiograph! que selon la revendication 1, caracterise en 
ce que la valeur d 1 attenuation attribuee a chaque 
volume est €gale a la somme des valeurs mesurees par 
ledit detecteur associe, divisee par le nombre des 

30 vignettes que contribuent S dormer ledit detecteur 
associe, et les valeurs d 1 attenuation des volumes sont 
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combinees par une combinaison numSrique sur des groupes 
(12) des volumes superposes aux differentes hauteurs de 
reconstruction. 

3) Proced6 de reconstruction d'une image 
radiographique selon la revendication l r caract^rise en 
ce que la valeur d' attenuation attribute a chaque 
volume est obtenue par retroprojection iterative des 
valeurs d 1 attenuation mesurees par les detecteurs (10), 
d es valeurs provisoires etant attributes aux volumes et 
corrigSes apr£s avoir ete projetees sur les detecteurs, 
en calculant des differences entre des sommes des 
valeurs provisoires sur des lignes de projection et les 
valeurs mesurees par les detecteurs sur lesdites lignes 
de projection, et en repartissant les differences sur 
lesdites lignes de projection pour corriger les valeurs 
provisoires . 

4) Proc€de de reconstruction d'une image 
radiographique selon la revendication 3, caracttrise en 
ce que les valeurs d 1 attenuation des volumes sont 
combinees en une somme sur des groupes (12) des volumes 
superposes aux differentes hauteurs de reconstruction. 
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